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1 Introduction 

Systematic approach to the multi-loop scattering amplitude study is now on the road, based 
on the fruitful progresses of tree and one-loop amplitude computations in the past a few 
years. After struggling with complicated calculations from Feynman diagrams for decades, 
the changes started with a new way of computing tree amplitudes. By complexifying 
the amplitude through certain momentum shifting on the complex plane, the amplitude 
becomes an analytic function of single complex variable A{z) with simple poles. The 
physical amplitude, defined as the amplitude at z = of the complex plane, is then 
calculated by the residues of A{z). The locus of residues are the phase space points where 
propagators are on-shell. Thus, the amplitude can be obtained by considering only on- 
shell diagrams with lower-point tree amplitudes as the input via Britto-Cachazo-Feng- 
Witten(BCFW) recursion relation [1, 2]. This point of view of scattering amplitude not 
only provides an efficient calculation method, but also greatly deepens our understanding 
of gauge field theories [6] . 

The breakthrough in tree amplitude calculations also inspires progresses in loop am- 
plitude calculations. Decades ago, the unitarity cut method [3] has already been used to 
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compute one-loop amplitude from tree amplitude [4, 5]. After the finding of BCFW re- 
cursion relation, this method is applicable to practical calculations [7, 8]. The difficulty is 
resolved because of the simple and compact tree amplitudes produced by the new method, 
so the focus of unitarity cut method is switched to the study of unitarity cut information. 

Take an one-loop amplitude as example, schematically it can be expanded on some 
one-loop integral basis [9], 



The basis Ij's are scalar integrals, i.e., integrals whose numerators of the integrand are 
1. They can be used universally for any one-loop amplitudes in renormalizable theories. 
Expansion coefficients q and the remainder 1Z are rational functions of external momenta. 
By cutting both sides of 1.1, i.e., computing the branch-cut discontinuities across various 
kinematical channels, we can get equalities between expressions in both sides, which are 
products of tree amplitudes. The coefficients can be determined by comparing expressions 
in both sides. Now the one-loop amplitude calculation is simple, because there are only 
finite number of integrals in the basis. For example, in 4-dimensional theory, the integral 
basis Ii can only be box, triangle, bubble or tadpole topologies, and the numerators are 
always constants. The solution space defined by the unitarity cut of these diagrams is 
simple: it is just the solution of one quadratic equation. In this sense, we can easily 
find a parametrization for the loop momentum, and work out the expansion coefficients Cj 
systematically from quadruple, triple, double and single cuts. Another way of extracting 
coefficients Cj is the generalized unitarity cut method [10, 11], which uses the information 
of cuts and also contour integration. 

The main concept of loop-amplitude calculations by unitarity-cut method is to expand 
the amplitude onto some known integral basis and then find the coefficients. So first of all, 
we need a set of integral basis. The basis should be large enough so that it is complete 
for expanding every loop amplitude, yet it should be as small as possible to simplify the 
calculation. It is somehow difficult to find the integral basis for multi-loop higher-point 
amplitudes. A conventional way is using Integrate-By-Parts (IBP) method [12] to find 
relations among different integrals. Sometimes the IBP calculation can be heavy, so instead 
we can define a set of integrand basis before integration. The denominator of the integrand 
is a product of propagators. The numerator is a polynomial of Lorentz Invariant Scalar 
Product (ISP). Ossola-Papadopoulos-Pittau(OPP) decomposition method [13] defines the 
integrand basis and the corresponding reduction scheme in the one-loop level. This method 
has been successfully applied to one- loop amplitude calculations, where the amplitude is 
decomposed into some Master Integrals (Mis) plus the rational terms [14-19]. 

It still takes some time to generalize above methods to multi-loop amplitude calcula- 
tions beyond the one-loop level calculation. The difficulty is obvious: the basis of multi-loop 
amplitudes is far more complicated than those of one-loop amplitudes. Some works have 
been done on the integral basis of diagrams such as two-loop double-box using IBP method 
[20] or generalized unitarity cut method [21-26]. By the latter method, the master con- 
tours, which produce the coefficients of double-box integrals in the basis decomposition, 
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are uniquely denned [23]. For more complicated diagrams, the study of integral basis is 
very complicated, thus we can consider the integrand basis first. After works in [27, 28], 
systematic study of integrand basis begins with the introduction of computational alge- 
braic geometry method, for example, Grobner basis, to multi-loop amplitude calculations 
[29, 30]. This method can determine a small enough yet sufficient set of the integrand basis, 
and also analyze the solution space of unitarity-cut equations systematically for two-loop 
[31-33] and also three-loop amplitudes [34]. 

The main objects we get, after applying the cuts on loop amplitudes, are unitarity- 
cut equations. These equations come from setting corresponding propagators on-shell, i.e., 
Di = lf — mf = 0, where i's ranges from all propagators being cut, li is the loop momentum 
and rrii is the mass of propagator. In this paper, we consider only massless theories, so 
rrii = 0. 

The studies on unitarity-cut equations are crucial for the application of unitarity meth- 
ods. For example, the unitarity-cut equations for AD massless four-point double-box dia- 
gram contain six irreducible components, and all the six components should be parame- 
terized to find the constraints for contours [21]. Furthermore, Simon and Kasper in [23] 
study the two- loop double-box diagram with six massive external legs, and find that the 
unitarity cut equations define an elliptic curve, i.e. a genus-one curve. This solution space 
can be parameterized by elliptic functions, but not rational functions. The structure of 
unitarity cuts for all other kinematics of 4D double-box can be easily illustrated via the 
degeneracy of the elliptic curve [23] . 

Unitarity cut equations may be simple, if the order of loop is low and the number of 
external legs is small. In these cases, one-variable complex analysis is sufficient for the 
analysis of unitarity cut equations. However, for high-loop diagrams or diagrams with 
many legs, the analysis would be complicated. Since the cut equations are all polynomial 
equations with complex coefficients, the natural tool to study them is algebraic geometry. 

The generalized unitarity cut of L-loop diagrams may define discrete points, complex 
curves or (hyper-)surfaces. All of them can be studied by algebraic geometry methods. 
In this paper, we consider the simplest non-trivial cases: the curves. More explicitly, we 
consider the cut equations of L-loop diagrams with 4L — 1 propagators in 4-dimensional 
theory. Since the topology of a complex curve is completely characterized by its genus, 
so we study the genera explicitly for the curves defined by generalized unitarity cut from 
one-loop to three-loop level, i.e. L = 1, 2, 3. The same mathematical approach should work 
for even higher-loop diagrams. However, the computation would be more involved. 

There are many ways to analyze genus and describe the global structure of multi-loop 
unitarity cuts, by algebraic geometry. In this paper, we mainly use the following methods, 

• The relation between arithmetic genus and geometric genus. In practice, we bira- 
tionally project the curve to a plane curve. For the latter one, the arithmetic genus 
can be easily found, and then we use the relation between arithmetic genus and 
geometric genus to find the genus of the original curve. 

• Riemann-Hurwitz formula. For a projection of a curve onto another one, this formula 
gives a relation between the genera of the two curves. 



-3- 



The genus of a curve is a birational invariant, i.e., invariant under rational re-parametrization 
Furthermore, if equations of unitarity cuts generate a genus-zero curve, then we can find a 
rational parametrization for the solution space. However, if the genus is larger than zero, 
the solution space cannot be described by rational parameters no matter how we choose the 
coordinates. So information of genus is the first judgement for the difficulty of an unitarity 
cut computation. 

Given a diagram, we first consider the kinematics with the maximal number of massive 
external legs. We call this case the prime case. Its unitarity cut sometimes gives high- 
genus curves. Then we consider the degenerate cases, i.e., the same diagram with fewer 
external legs or with several massless external legs. Then usually we get a reducible curve, 
which is the union of several irreducible curves. Each curve's genus is lower than the one 
of the prime case. For the double-box case in [23], we show again that the global structure 
of unitarity cut is given by the prime case, i.e., if we find the intersection of curves from 
the degenerate case and "sew" them together, we get the topology of the prime case. In 
this sense, the genus information is the guideline for the global structure of unitarity cut 
solutions for all degenerate kinematics. 

Throughout this paper, we find that the genus of high-loop unitarity cuts can be higher 
than 1. For example, for a generic non-planar two-loop crossed-box diagram with 6 massive 
momenta, the on-shell equations of unitarity cuts generate a complex curve with genus 3. 

The paper is organized as follows. In section 2, we present some mathematical back- 
ground, where basic concepts of algebraic curves, such as genus and singular points, are 
introduced. In section 3,4 and 5, we study the genera of algebraic curves from one, two 
and three-loop diagrams respectively. Detailed analysis is given for two-loop diagrams, and 
some typical diagrams of three-loop are also carefully studied. Discussions and comments 
are presented in the conclusion. 

2 Mathematical preliminaries 

In this section, we review the definition and computation of the genus of a complex algebraic 
curve. For detailed mathematical proofs, please refer to [35] [36]. 

There are two types of genera for complex curves, arithmetic genus and geometric 
genus. The topology of a complex curve is characterized by its geometric genus. Roughly 
speaking, the geometric genus is the number of handles of a complex curve. (The rigorous 
definition of the two genera will be given in this section.) Throughout the paper, if not 
specified, "genus" means geometric genus. 

We calculate the geometric genus of a complex curve directly by computational alge- 
braic geometry. The strategies are, 

• Calculate the arithmetic genus of a complex curve. This is easily done by polyno- 
mial computations, which is automated by the computational algebraic geometry 
softwares, for example, 'Macaulay2' [37] . If the curve is non-singular, then the arith- 
metic genus equals the geometric genus and we are done. Otherwise, we project 
the curve onto to a plane curve and re-calculate the arithmetic genus. For the plane 
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curve, the geometric genus is simply related to the arithmetic genus by singular-point 
counting and the blow-up process. 

• Alternatively, project the curve onto a known curve with smaller genus. This projec- 
tion needs not to be birational. Then we use Riemann-Hurwitz formula to determine 
the genus of the first curve. 

2.1 Projective algebraic curves and the arithmetic genus 

Consider a complex curve C in C n , defined by polynomial equations 

fi(z 1 ,...z n ) = ... = fk(zi, ...z n )=0. (2.1) 

The polynomials fi,...fk generate an ideal J = • • • fk)- In the formal language, C is 
the zero locus of J. 

To discuss its topology, we need to extend C to a projective curve C in CP n . Define 
homogenous coordinates by Z{ = Zi/Zo, C is defined by homogenous polynomials 

Z^ f -f&,...^)=F(Z Q ,...,Z n ) = 0, V/6J. (2.2) 

All such P's generate a homogenous ideal I. In practice, to get /, we do not need to 
find the homogenous form of all polynomials in J. Instead, it is sufficient to consider the 
homogenous form of polynomials in the Grobner basis of J, in a graded monomial order 
[38]. 

Furthermore, we assume that C is irreducible, i.e., cannot be written as the union of 
two different curves. Equivalently, it means that I is a prime ideal [35]. For the cases 
of reducible curves, we consider the genus of each irreducible component. Irreducible 
components can be found by primary decomposition [35], 

I = f]li, i = l,...,k, (2.3) 

i 

where all ij's are prime ideals. Primary decomposition technique has been used in high- loop 
unitarity cuts [29] [34] [31]. 

The arithmetic genus can be defined by the Hilbert polynomial of the quotient ring 
S/I, where S = C[Zq, Z\, . . . Z^\. S/I is a graded S-module, 



5// = 05i. (2.4) 



i=0 

There exists a unique polynomial P(x) and a positive integer N such that, for n > N, 
n E Z, 

P(n) = dime S n . (2.5) 

This polynomial is called Hilbert polynomial (of the projective curve C). Then the arith- 
metic genus of C is defined to be 

g A {C) = -P(0) + 1 . (2.6) 
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In practice, the Hilbert polynomial and arithmetic genus can be easily calculated by 
Grobner basis method. In particular, for a plane curve, i.e. a curve in CP 2 , defined by one 
homogenous polynomial with degree d, the arithmetic genus is simply [35] 

gA = l(d-l)(d-2) . (2.7) 
2.2 Singular points and the geometric genus 

The arithmetic genus is not directly related to the topological properties of a curve, since 
intuitively it counts not only the number of handles but also the singular points. A singular 
point on C is a point (ao,oi, ....a n ) such that the rank of the Jacobian matrix, 



dFi. 

(a , .. 



1 < i < k , < j < n , (2. 



9G = 9A~ 

PeSing(C) 



dZj 

is less than n — 1 [35], where Fi, 1 < % < k are the generators of equations for the curve. A 
singular point on C is normal if all tangent lines on the singular point are distinct. 

For a smooth curve, i.e., an irreducible projective curve without a singular point, we 
define the geometric genus to be its arithmetic genus [35]. For an irreducible projective 
curve C with singular points, there exists a smooth irreducible projective curve C, which is 
the normalization of C [35], 

■k:C^C . (2.9) 

We define the geometric genus of C to be the arithmetic genus of C. This definition is 
consistent with the topological definition of the genus. In particular, for an irreducible plane 
curve with normal singular points only, the geometric genus is related to the arithmetic 
genus as 

\Mvp - 1) (2-10) 

Sing(C) 

= I(d-l)(d-2)- ]T i M p(A*P-l), (2-11) 

PeSing(C) 

where Sing(C) is the collection of all singular points P on C. \xp is the multiplicity of P, 
i.e., the number of tangent lines at P. In practice, Sing(C) can be found by computational 
algebraic geometry softwares, like Macaulay2 [37]. 

We may frequently re-parametrize a curve, and intuitively, the topological information 
should be invariant after re-parametrization. To realize the re-parametrization process rig- 
orously, we introduce the concept of birational map: A rational map r from the irreducible 
curves, C\ and C2, is a map from a non-empty open set of C\ to C2, such that, in terms of 
the coordinates, all components of / are rational functions. A birational map / between 
two irreducible curves is a rational map / from C\ to C2 with an inverse rational map g 
from Ci to C\, such that gf is the identity map on an open set of C\ and fg is the identity 
map on an open set of C2 ■ The geometric genus is a birational invariant [35] . 

In particular, since CP 1 has the genus zero, a project curve with go > is not birational 
to CP 1 . It implies that there is no rational parametrization for such a curve. 

We present some simple examples on the genera of algebraic curves: 
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1. y 2 = x 3 + 1. This is a cubic plane curve. The corresponding projective curve is 
Y 2 Z = X 3 + Z 3 . There is no singular point on the projective curve. The arithmetic 
genus is (3 — 2) (3 — l)/2 = 1 and the geometric genus is also 1. So this is an elliptic 
curve. 

2. y 2 = x 3 + x 2 . This is also a cubic plane curve. The corresponding projective curve 
is Y 2 Z = X 3 + X 2 Z. There is a singular point at X = 0, Y = 0, Z = 1 and there 
are two distinct tangent lines at this point: x — y = and x + y = 0. So this is an 
ordinary singular point with the multiplicity 2. The arithmetic genus is still 1 but 
the geometric genus is 1 — 2(2 — l)/2 = 0. So this is not an elliptic curve. We can 
blow up the singular point at (x,y) = (0,0) as follows [35]: define a new variables 
t via y = xt, and the equation becomes — x 2 (l — t 2 + x) = 0. Factorize it and keep 
the curve part, we get a new curve 

1 - t 2 + x = 0. This new curve is birational to 
the original curve, and is smooth. Since the new curve is conic, so we show that the 
original curve is birationally equivalent to a conic and can be rational parameterized. 

3. The projective curve, 

XT-YZ = 0, X 2 Z + y 3 + YZT = 0, XZ 2 + Y 2 T + ZT 2 = (2.12) 

is non-singular [36] . From Macaulay2, we find that the arithmetic genus is 2, so its 
geometric genus is also 2. 

There are many different ways to calculate the genus of an algebraic curve. In this 
paper, for the algebraic curves from unitarity cuts in two-loop and three-loop orders, we 
mainly use the relation between arithmetic genus and geometric genus for the computation. 
Explicitly, if the curve is non-singular, then we simply calculate its arithmetic genus and 
by the definition, it is the geometric genus. Otherwise, we birationally project the curve 
onto a plane, and then use (2.10) to determine the genus if all singular points are normal. 
If some singular points are not normal, we use the blow up process to resolve them, as 
described in [36]. 

Anther useful formula for curve genus computation is Riemann-Hurwitz formula [35]. 
Let / : C\ — > C2 be a covering map between two irreducible complex curves C\ and C2, 
and n is the degree of /. Then g{C\) and g(C2), the genera of C\ and C2 are related by 

2g(d) - 2 = n ■ (2g(C 2 ) - 2)) + £ (e P - 1) , (2.13) 

PeCi 

where ep is the ramification index of P. For all but finite points on C%, the ramification 
index equals one so it is just a finite sum. Therefore if g{C2) is known, we can get g(C\) 
by counting the ramified points on C\. 

3 Genera of algebraic curves from one-loop diagrams 

For one-loop diagrams, there is only one loop momentum, and it has 4 components. Di- 
agrams with 3 propagators after maximal unitarity cut define algebraic curves. The only 
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diagram we should consider is the triangle diagram. The propagators are given by 

D = £ 2 , D 1 = (l-p 1 ) 2 , D 2 = {l- Pl -p 2 f , (3.1) 

where external momenta at each vertices are pi,p 2 ,P3, an d we assume they are all massive. 
We can get 2 polynomials L\ = Dq — D\, L 2 = Dq — D 2 , which are linear in I. If we use 
4 variables (x±, x 2 , Xs, X4) to expand the loop momentum, then L\ = L 2 = are 2 linear 
equations in them. By solving these equations and substituting solutions back to Dq, we 
get a quadratic polynomials of 2 variables. The resulting equation Dq = is then a conic 
section, and topologically equivalent to a genus-0 Riemann sphere. 

If some of the external momenta are massless, the on-shell equations become degener- 
ate. The topological picture can be given by combing another Riemann sphere from the 
original one, and the result is that two Riemann spheres are connected at a single point. 
This agrees with the analysis of intersection pattern in [31]. 

4 Genera of algebraic curves from two-loop diagrams 

In this section, we calculate the genera of all curves from the 4D two-loop unitarity cuts, 
for massless theories 

For 4-dimensional two-loop diagrams, there are 8 components from two loop momenta. 
The 8 variables associated to the 8 components of loop momenta could be defined as the 
expansion coefficients of loop momenta on a chosen momentum basis. They can also be 
defined by Lorentz invariant scalar product of loop momentum and external momentum. 
The resulting propagators, which expressed as polynomials of these 8 variables, will be 
different according to the definitions of momentum basis and variables. But in general, 
they are always quadratic polynomials and only differ by a linear transformation. 

If we apply the maximal unitarity cut on the diagram, all propagators become on-shell, 
then we get a polynomial equation system. These equations define an algebraic set, i.e. 
the solution set of cut equations in the space spanned by 8 variables. When considering 
diagrams with 7 cuts, the solution set is described by 7 equations with respect to 8 variables. 
It is an algebraic curve. 

It is not easy to analyze this curve directly from its defining polynomial equation 
system. So the first step is to eliminate as many variables as possible from the linear 
equations which can be constructed from 7 on-shell equations. In practice, it is always 
possible to obtain 4 linear equations. By solving them, we can express 4 variables as linear 
combinations of the remaining 4 variables. The resulting 3 quadratic polynomial equations 
of 4 variables is isomorphic to the original algebraic set. Furthermore, by coordinate 
transformations, we can rewrite the 3 equations as a single meromorphic function of 2 
variables. It defines a plane curve, which is birationally equivalent to the original algebraic 
set. Though the explicit definitions of curves in above three descriptions are different, the 
topological invariant objects remain the same. So we can get a unique result for the genus. 
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Figure 1. Two- loop diagrams with 7 propagators: (a) planar pentagon-triangle diagram, (b) planar 
double-box diagram, (c) non-planar crossed-box diagram. All external momenta are out-going and 
massive. The loop momenta are denoted by £\^£i. 

4.1 Two-loop planar pentagon-triangle diagram, genus 

As a warm-up, we first analyze this simple diagram, in Figure. (1. a). The 7 propagators 
are given by 

D = ll, D 1 = (£ 1 - Pl ) 2 , D 2 = {£ l -p l -p 2 ) 2 , D 3 = (£ 1 -p 1 -p 2 -p 3 ) 2 , 
D = £ 2 2 , L>i = (£ 2 — P4) 2 , Do = (h+£ 2 +p 5 ) 2 . (4.1) 

There are 4 propagators containing £\ only, and from the on-shell equations D{ = 0, we 
can completely fix the variables in £\. These equations Di = 0, i = 1,2,3,4 have two 
solutions. For each solution, since £\ has already been completely fixed, we can treat it as 
a four-vector similar to an external momentum. Thus the on-shell propagators highlighted 
by dashed lines in Figure (l.a) can be effectively treated as massless external momenta. 
Then for each solution, the result of two-loop planar pentagon-triangle diagram is the same 
as the result of the one-loop triangle diagram. It is associated with a genus-0 Riemann 
sphere. In some specific kinematic configurations where degeneracy exists, the topological 
picture becomes that of two Riemann spheres connected at a single point. 

4.2 Two-loop planar double-box diagram, genus up to 1 

It is already known in [23] that on-shell equations of the two-loop double-box diagram, 
with six massive legs, define an one-dimensional elliptic curve which is associated with a 
genus- 1 torus. 

We can reproduce this result by directly computing the arithmetic genus and singular 
points of the curve. Taking the convention shown in Figure (l.b), the 7 propagators can 
be written as, 

D = £ 2 , D 1 = {£ 1 -p 1 ) 2 , D 2 = (£ l -p l -p 2 ) 2 , 

D = £\ , D 1 = (£ 2 - pa) 2 , D 2 = (£ 2 -pa- p A f , D = (h +£ 2 + p 5 ) 2 . (4.2) 

In order to compute the arithmetic genus and singular points, we first rewrite above 
equations by expanding all momenta on a momentum basis. A simple analysis is as follows. 
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I 

(a) (b) (c) (d) 

Figure 2. Topological pictures of on-shcll equations from the two- loop double-box diagram under 
specific kinematic configurations. The pictures should be understood as complex curves, or two- 
dimensional real surfaces, (a) the curve is irreducible and the solution set is a torus (b) the curve 
has 2 irreducible branches, (c) the curve has 4 irreducible branches, (d) the curve has 6 irreducible 
branches. For general kinematics the curve is genus 1. In degenerate limit, tubes shrink to points 
along dashed lines. The resulting Ricmann surfaces for each branch can only be a sphere. This 
explains why we get Riemann spheres connected by points and linked in a chain in degenerate 
limits. 



We have variables (xi, x 2 , xs, x±) for 4 components of l\ and {yx,yi, 2/3, 2/4) for 4 components 
of £2- Using on-shell equations Di = 0, Di = 0, we can eliminate 2 variables of Xi and 2 of yi, 
and get 3 quadratic polynomials Q\ = Qi(xi,x 2 ), Q2 = (32(2/1,2/2), Q3 = Qsi^i, x 2 , yi, 2/2)- 
This can be done by Yang's BasisDet package[29]. We can further eliminate 2 variables 
and 2 equations via Grobner basis method and get a plane curve. This plane curve has 
degree 8, so the arithmetic genus is 

(d-l)(d-2) 

9A = g = 21 • ( 43 ) 

There are 10 singular points, of which 8 have the multiplicity fip = 2 and 2 have the 
multiplicity \ip = 4, so the geometric genus is 

9G = 9A~ YI 2 flp ( fip - 1 ) = 21 - 8xl - 2x6=1 - ( 4 - 4 ) 

PeSing(C) 

This result is consistent with that in [23]. 

Another explicit way of getting an equivalent plane curve can be taken as follows. 
Through coordinate transformation it is easy to rewrite Q\ — > = x^x^ — c±, Q2 — > 
Q'2 = 2/12/2 — c 2> so we can do the following substitution x[ = ci/x 2 ,y[ =£2/2/2 i n ^3- The 
resulting Q' s = n(x' 2 , y^jdix^ 2/2) is a meromorphic function which is equivalent to original 
equations, and the numerator n(x 2 ,2/2) defines a plane curve. This plane curve n(x' 2 ,2/2) 
has the degree 4, so we have gA = 3. There are 2 singular points of the multiplicity /ip = 2, 
thus we get the geometric genus go = 3 — 2 = 1. This result again agrees with that given 
in [23]. 

The topological pictures of genus of degenerate algebraic curves under specific kine- 
matic configurations can be easily found, because in degenerate limit the tube shrinks to 
point, and there is only one handle for the torus. The contraction will break the handle, and 
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we should get genus-0 Riemann spheres connected at points. Equations of each irreducible 
branch can be obtained by Macaulay2 [37] via primary decomposition. By computing the 
arithmetic genus and singular points, we found that their geometric genera are indeed zero. 
The 3 kinds of contraction shown in Figure (2) describe topological pictures of degenerate 
algebraic curves under all possible specific kinematic configurations. There are 2, 4 or 6 
Riemann spheres connected by points along dashed lines where tubes shrink, and linked 
adjacently into a chain. This is exactly the picture described in [23]. 

4.3 Two-loop non-planar crossed-box diagram, genus up to 3 

For the two-loop non-planar crossed-box diagram, the explicit form of propagators for 
generic kinematics is very complicated. Abstractly, we write the 7 propagators as, 

D = £j, £>i = (4 -Pi) 3 , D 2 = (£ l - Pl - P2 f , D =£ 2 , 

Di = (£2 - Psf , D = (h+£ 2 + P5 ) 2 , D 1 = (£ 1 +£ 2 + P4 + P5 ) 2 . (4.5) 

From above propagators, we can at most get 4 polynomials which are linear in £, 

L\ = Dq — D\ = 2£ lPl -pj , L 2 = D -D 2 = 2£ 1 (p 1 + P2 ) - { Pl + P2 f , 
L 3 = D - D x = 2£ 2P3 -p 2 , L A = D -D 1 = -2{h + £ 2 + P5 )p 4 - P \ , 

and 3 polynomials which are quadratic in £, 

Q 1 = D = £ 2 , Q 2 = D =£ 2 , Q 3 = D - D - Do = 2h£ 2 + 2(h + £ 2 ) P5 + P 2 . 

If we use momentum basis ei,e 2 ,e%,e4 to expand loop momenta as, 

£1 = x 2 ei + xie 2 + x 4 e 3 + x 3 e 4 , £ 2 = y 2 ei + y±e 2 + y 4 e 3 + y 3 e 4 , (4.6) 

and the 8 variables are taken as the expansion coefficients (x±, x 2 , X3, x 4 ) and (2/1,2/2,2/3, 2/4), 
then there will be quadratic terms from products of two loop momenta £{£j, and linear 
terms from products like l% P j. Systematically we can generate the Grobner basis by 7 
on-shell equations Di = D, t = D, t = 0, and find all linear terms in the Grobner basis. 
We can remove as many variables as possible by solving these linear terms, and express 
other terms with remaining variables. In this case, we get an equivalent but simplified 
polynomial equation system. A more explicit way is to solve linear equations Li = 0, and 
express 4 variables as linear functions of remaining 4 variables. Substituting solutions back 
to quadratic polynomial equations Qi = 0, we get an alternate algebraic set. It is possible 
to analyze the three quadratic equations by algebraic geometry program such as Macaulay2 
[37]. 

Let us consider the two-loop non-planar crossed-box diagram drawn in Figure (l.c). All 
external momenta are general, and we construct momentum basis &i from pi,P3 . External 

1 Momentum basis a have the following orthogonal and normalization relations eie2 = 6364 = 1, other- 
wise 0. The explicit form can be found by Gram-Schmidt process. 
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momenta are expanded as 

4 

pi = a\e\ + Qi 2 e 2 , pi + P2 = ^ a2iei > ^ = ^ liei + ^ 12e 2 > 

i=l 

4 4 

P5 = ^ Tiiei , P4 + f>5 = 72iei , Pe = -(Pi + P2 + Vi + Pi + Ps) , (4.7) 
1=1 1=1 

where a, /?, 7 are the projection coefficients of momenta on the corresponding momentum 
basis, for example an = pie2, ai2 = pie±. Using this expansion, the four linear polynomials 
Li become 

L\ = 2(anxi + a 12 X2 — ana^) , 

L 2 = 2(a 2 ixi + a 22 x 2 + a 2 3X 3 + a 2 4^4 — a 2 ia 22 — a 23 a 2 4) , 

L 3 = 2(p uyi + £122/2 - faiPvt) , 
4 

La = -2 ^(xi + yi)(rf2i - 7it) - 2(721722 + 723724 - 7ii7i2 - 7137m) , (4.8) 
i=i 

and the three quadratic polynomials Qi become 

Qi = 2(xix 2 + X3X4) , Q 2 = 2(yiy 2 + y 3 y 4 ) , 

4 

Q 3 = 2(xiy 2 + x 2 yi + x 3 y 4 + x 4 y 3 ) + 2 ^J(x* + yi)7i» + 2(711712 + 713714) • (4.9) 

i=l 

The equations until now are still simple, but after solving Lj = and substituting solu- 
tions back to Qi, they will become quadratic polynomials with complicate coefficients of 
a,/3,7. Anyway, Qi's are still 3 quadratic polynomials in 4 variables, and they define a 
one-dimensional curve. 

As we have mentioned, different definitions of momentum basis and variables, and 
the choices of choosing remaining variables after solving linear equations, will describe 
different curves. But they are all birationally equivalent. Birational invariants will stay 
the same in each description. So we project the three quadratic equations onto a plane 
curve. By a coordinate transformation, we can eliminate 2 equations as well as 2 variables, 
and the resulting equation defines a plane curve describing an equivalent algebraic set. 
This equation has the degree d = 8, so the arithmetic genus is, 

g A = I(d-l)(d-2) = 21 . (4.10) 

The singular points can be obtained directly from the definition, i.e., for projective plane 
curve defined by homogeneous polynomial P(x, y, z) of degree d, the singular points are 
solutions of equations 

P(x,y,z) = P^x,y,z) = P'(x,y,z) = P' z (x,y,z) = . (4.11) 
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For this plane curve, there are 18 normal singular points of the multiplicity fip = 2, so the 
geometric genus is again 



Another way of getting a plane curve from original on-shell equations is to use the 
elimination process via Grobner basis method. We can use BasisDet package to generate 
on-shell equations, where the variables are defined as independent Lorentz Invariant Scalar 
Products (ISPs) of loop momentum and external momentum. Using these cut equations 
and Grobner basis associated with them, we get a plane curve. The degree of the plane 
curve is 8, so we have gA = 21. There are 18 normal singular points of the multiplicity 2, 
so again we get the geometric genus go = 21 — 18 = 3. 

To summarize, the maximal unitarity cut of two-loop non-planar crossed-box diagram 
can fix 7 components of loop momenta, and the solution set of on-shell equations is de- 
scribed by a free parameter corresponding to the remaining one degree of freedom. The 
variety is associated with a genus-3 Riemann surface. 

4.3.1 Degeneracy under specific kinematics 

In the prime case, external momenta pi are massive and the algebraic set defined by on-shell 
equations is irreducible. Under some specific kinematic configurations, the variety will be 
decomposed to many irreducible ideals after primary decomposition. For example, if p\ is 
massless in Figure (l.c), the algebraic set decomposes into 2 irreducible ideals, i.e., it has 
two branches. Each irreducible branch is associated with a Riemann surface, and it should 
have connection with the original genus-3 Riemann surface. 

We can compute the geometric genus for each irreducible branch using the same method 
in previous section. Let us start from the kinematic configuration where at least one 
momentum of p±,P2 is massless. Simple calculation shows that for pi = 0, we have Q\ = 
so Q\ = implies that either X3 = or X4 = 0. For p\ = 0, we have Q\ = 
fi(x3, X4)f2(x3, X4), where /i,/2 are linear polynomials of X3,x^, and Q\ = implies that 
either /1 = or /2 = 0. Combined with Q2 = 0, Q3 = 0, there are two branches. For each 
branch, we can transfer the 3 equations to a plane curve by coordinate transformation. 
This plane curve has the degree 4, so the arithmetic genus is gA = 3. There are 2 normal 
singular points of the multiplicity /ip = 2. Finally the geometric genus is go = 3 — 2 = 1 for 
each branch, so they are associated with a genus-1 torus. The two tori intersect at 2 points 
[31]. From these it is easy to obtain that, the topological picture under this kinematic 
configuration is given by contracting two tubes to points along dashed lines, as shown in 
Figure (3.b). The up and down parts become tori, and they are connected by 2 points, 
which are originally 2 tubes. 

If we consider kinematic configurations where only one of P3,P4 is massless, or only 
one of P5,pe is absent (which means that only one vertex of those attached to P3,P4,P5,Pe 
is a three- vertex with chirality), the algebraic set is reducible, and it has two irreducible 
branches. For p\ = 0, we have Q2 = J/3J/4, and we can get two branches directly from 
Q2 = 0. For other kinematic configurations, it is not easy to see the degeneracy directly 
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Figure 3. Topological pictures of degenerate on-shell equations from non- planar two- loop crossing- 
box diagram under specific kinematic configurations, (a) the curve is irreducible, the geometric 
genus is 3, (b) the curve has 2 irreducible branches, each branch is genus 1, and they are connected 
by two points, (c) the curve has 2 irreducible branches, each branch is genus 0, and they are 
connected by four points, (d) no kinematic configuration corresponds to this topological picture. 

from on-shell equations, since none of Qi,Q2,Q3 can be obviously factorized into two 
parts. Then we generate an ideal from Q±, Q2, Q3, get two irreducible ideals using primary 
decomposition method of Macaulay2 [37]. For each irreducible ideal, we can obtain a 
plane curve after coordinate transformations. Each plane curve has the degree 4, so the 
arithmetic genus is gA = 3. There are 3 normal singular points of the multiplicity [ip = 2, 
so we have gc = 3 — 3 = 0. It means that for all kinematic configurations considered 
in this paragraph, each irreducible branch is associated with a genus-0 Riemann sphere. 
Another way of analysis is to get a plane curve directly from Qi = Q2 = Q3 = 0, and the 
resulting equation has the degree 8. This equation can be factorized to 2 polynomials of 
the degree 4, and each factor is equivalent to one irreducible branch. We can again get 
geometric genus go = 3 — 3 = 0. These two Riemann spheres intersect at 4 points [31], so 
the possible topological picture can be given by contracting 4 tubes to points along dashed 
lines as shown in Figure (3.c). The left and right parts become Riemann spheres, and they 
intersect at 4 points. 

Above discussion involves all kinematic configurations where the algebraic set has 
two irreducible branches. From a genus-3 Riemann surface, there is another possibility 
of contracting two tubes to get two Riemann surfaces as shown in Figure (3.d), and the 
resulting topological picture is: a genus-0 Riemann sphere and genus-2 Riemann surface 
intersect at two single points. But no kinematic configuration is found to guarantee the 
on-shell equations having such degeneracy. Intuitively, they can be understood as the 
consequence of symmetry between two branches. Usually, the two branches are related by 
the parity symmetry and they would have the same genus. 

If we denote kinematic condition as, 

• K\\ at least one of pi,P2 is massless, 

• K% : t?3 is massless, 
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• K\: is massless, 

• iff : P5 is absent, 

• K\: pq is absent, 

then the general rules for degenerate topological picture can be given as follows: condition 
K\ will contract the tubes at dashed lines in Figure (3.b), and condition K%, K\, K% or 
K\ will contract the tubes at dashed lines in Figure (3.c). Topological picture of cases with 
more than two irreducible branches can be determined by combining these conditions. 

Explicitly, if we combine one condition K\ with one of , K\, K$, K^, the topolog- 
ical picture is given in Figure (4. a), which is the overlap of (3.b) and (3.c). This gives 4 
Riemann spheres, which means that there are 4 irreducible branches. The on-shell equa- 
tions of this kinematic configuration indeed have 4 irreducible branches, and by the explicit 
computation, each one is a genus-0 curve, which agrees with the topological picture. 

If we combine one condition of K^K^ with one of K^,K^, the topological picture is 
given in Figure (4.b). This gives 4 Riemann spheres. It is not hard to find that Figure (4.c) 
is given by further combining K\ with above configuration, and the resulting picture gives 
6 Riemann spheres. For kinematic configurations with both K^^K^ or both K^,K^, the 
picture is given by overlapping the double copy of Figure (3.c), as shown in Figure (4.d). 
This also gives 6 Riemann spheres. 

Furthermore, combined with K\, we get the topological picture (4.e), which has 8 
Riemann spheres. (This case was studied in [28].) We computed the geometric genera of 
all irreducible branches in this paragraph, and found that they are all zero. Intersection 
points computation from on-shell equations of irreducible branches, also agrees with this 
analysis. This verifies the result in [31]. 

5 Genera of algebraic curves from three-loop diagrams 

We also consider algebraic curves defined by equations of maximal unitarity cut from 4- 
dimensional three-loop diagrams. There are 4 x 3 = 12 components in loop momenta, so 
diagrams with 11 propagators would generate algebraic curves. We study the genera of 
several three-loop curve examples. Although the algebraic system is more complicated in 
these cases, the genus computation process is similar to that of two-loop cases. We will 
project curves to plane curves, and use the knowledge of arithmetic genus and singular 
points to compute the geometric genera. 

5.1 Three-loop planar pentagon-box-box diagram, genus up to 1 

We consider the generic pentagon-box-box diagram with 9 massive external legs, as shown 
in Figure (5. a). The 11 propagators are given by 

D =£j, D 1 = (£ 1 -p 1 f , D 2 = (h - Pl - p 2 f , D 3 = (4 - Pl - P2 - p 3 f , 
D = l\ , D 1 = (£ 2 -p 5 ) 2 , D 2 = (£ 2 - P4 - P5 f , D = t\ , 

D 1 = (£ 3 - Pl -p 2 -p 3 -p 6 -p 7 ) 2 , Do = (4 - h - pg) 2 , D 1 = (e 3 + e 2 + P8 ) 2 . 
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(d) (e) 

Figure 4. Topological pictures of degenerate on-shell equations from non-planar two-loop crossing- 
box diagram under specific kinematic configurations, where the solution set has more than 2 
branches. These 5 pictures includes the degeneracies under all possible kinematic configurations. 
Each irreducible branch is genus sphere, and they are connected by points along dashed lines 
where tubes have been contracted. 



P3 



P2 




Figure 5. Planar three-loop diagrams with f f propagators: (a) pentagon-box-box diagram, (b) 
box-pentagon-box diagram. All external momenta are out-going and massive. The loop momenta 
are denoted by t\^2,(-z- 

The 11 on-shell equations Di = Di = Di = Di = of maximal unitarity cut define 
an algebraic curve, and its topology is actually simple. To see this, note that 4 on-shell 
equations Di = contain only £i, and we can solve them first. Since t\ has 4 components, 
these 4 equations will completely fix l\. Explicitly, after solving 3 linear equations from 
Di = 0, the remaining one becomes a quadratic equation of a single variable. This equation 
will give two solutions for £\, and on each solution, i\ is completed fixed. Hence i\ can be 
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treated as a 4- vector similar to external momenta. The remaining 11 — 4 = 7 cut equations 
for l 2 and l 3 are the same as the two-loop double-box case. 

More explicitly, for each solution of £%, the on-shell propagators highlighted by dashed 
lines in Figure (5. a) can be treated as massless external momenta. So the discussion reduces 
to the two-loop double-box diagram with general kinematics if pe,P7 are not absent. As we 
know it is a genus-1 torus. In some specific kinematic configurations of (pi,P5,P6,P7,P8,P9), 
the algebraic curves are degenerate, and the topological pictures are given by genus-0 
Riemann spheres connected at points and linked in a chain, just as the two-loop double- 
box diagram. The two solutions of £\ are separated. 

5.2 Three-loop planar box-pentagon-box diagram, genus up to 5 

On-shell equations of the three-loop box-pentagon-box diagram are much more compli- 
cated. Naive elimination of linear equations will generate 5 quadratic equations in 6 vari- 
ables, and they define a one-dimensional curve. More explicitly, we consider the generic 
box-pentagon-box diagram with 9 massive external momenta as shown in Figure (5.b). 
Among the 11 propagators, there are 3 containing l\ only 

D =£l, £>! = (h - pif , D 2 = {£ x - Pl - p 2 ) 2 , (5.1) 
3 containing £ 2 only 

D = l\ , D x = (£ 2 - p 5 ) 2 , D 2 = (£ 2 - p 4 - P5) 2 , (5.2) 

and 3 containing £ 3 only 

Do=l\ , Di = (£ 3 -pi-p 2 -p e - p 7 f , D 2 = (£ 3 -pi-p 2 -p3-p e - p 7 ) 2 ■ (5.3) 

The remaining two propagators contain terms of mixed loop momenta, and they are given 
by, 

A, = (£ 3 ~h- Psf , D x = (£ 2 +£ 3 + p 8 ) 2 . (5.4) 

We have 4 variables for each loop momentum, denoted as (x±, x 2 , x 3 , X4) for £±, (yi,y 2 ,y 3 , 7/4) 
for £ 2 and (z±, z 2 , z 3 , z^) for £ 3 . As usual, for each loop momentum we have 2 linear equa- 
tions, and by solving them, we can eliminate 2 variables and 2 equations. So finally we get 
5 quadratic polynomials, 

Qi(x 3 ,x4) = D , Q 2 (y 3 ,yi) = Do , ^3(^3, u) = Do , 

Qi(x 3 ,x 4 , 03, za) = Do - D - D , Q 5 {y 3 , y 4 , z 3 , z 4 ) = D ± - D - D . (5.5) 

Notice that Qi(x 3 , X4) = 0, Q 2 {y 3 -,yi) = and Q 3 (z 3 , £4) = are conies. So it is always 
possible to find rational parametrization for each loop momentum, and express them by 
one free parameter as l\ = £\{x), t 2 = £ 2 (y) and £ 3 = £ 3 (z) 2 . The remaining two equations 

2 By coordinate transformation, we can rewrite a non-degenerate conic Q(x,y) = in the form 
Q'(x',y') = x'y' — c = 0, and express one parameter as a rational function of the other as x' — c/y'. 
So the loop momentum can be expanded as a rational function of a single variable. 
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Dq and D\ are meromorphic functions after substituting li back, and the numerators define 
an algebraic plane curve as, 

a 2 (x)z 2 + a\(x)z + a (x) = , 
b 2 (y)z 2 +b 1 (y)z + b (y) = . (5.6) 

For generic massive external momenta, and bi are quadratic polynomials. We can elim- 
inate z by computing the resultant of these two equations to get a plane curve F(x, y)=0- 
This plane curve is birationally equivalent to (5.6) via the inverse map, 

z _ -a 2 {x)b {y) + ao(x)b 2 {y) , g _ 

a 2 (x)b 1 (y) - ai(x)b 2 (y) 

so they have the same geometric genus. For generic kinematic configurations, the plane 
curve F(x,y) has the degree 8, so the arithmetic genus is 

^- P-'K'-q -H. (5.8) 

There are 4 singular points of the multiplicity \ip = 2, and 2 singular points of the multi- 
plicity /ip = 4, so the geometric genus is given by 

9G = 9A~ \m»p-V = 21 -4x1-2x6 = 5. (5.9) 

PeSing(P) 

So this is a genus 5 curve. 

Alternatively, we can also find the genus of (5.6) by Riemann-Hurwitz formula (2.13) 3 . 
Let C be the complex curve defined by (5.6). The projection / : (x,y,z) H> z is a covering 
map from C to the complex plane of z. For all but finite points, this is a four-fold covering 
so the degree of the map is 4. We determine the ramified points as follows: 

1. Consider the first equation in (5.6) as a quadratic equation for x and the discriminant 
is A\(z). For a fixed value of z, there are two corresponding x unless A\(z) vanishes. 
We can checked that Ai(z) = have four distinct roots on the z plane, namely z±, 
z 2 , Z3 and z^. 

2. Consider the second equation in (5.6) as a quadratic equation for y and the discrim- 
inant is A 2 (z). We can checked that A 2 (z) = have four distinct roots on the z 
plane, namely z$, zq, zj and z%. 



3. We checked z^s are all disctint, i = 1, . . . , 8. Hence for each Zi, i = 1, ... 4, there 
are two corresponding points on C, namely (xi,y\ , zi) and (x« , y • , Z\ ) . For each 
Zi, i = 5,... 8, there are two corresponding points on C, namely (x^\ yi, Z{) and 
(x- ,yi, Zi). So there are 2x4 + 2x4 = 16 ramified points on C and each point has 
the ramification index 2. 



3 We thank Simon Caron-Huot for the discussion of using this formula. 
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Then we can use Riemann-Hurwitz formula (2.13) and the fact that the complex plane has 
the genus zero, 

25(C) - 2 = 4 x (2 x - 2) + 16 • (2 - 1) , (5.10) 

Again we get the same conclusion that the curve has the genus 5. 

Finally, we can also find the genus with the help of two- loop-diagram genus information. 
If we forget the cut equations involves I2 for a while, the rest equations form the hepta- 
cut equations of the two-loop massive double box diagram. Intuitively, it is clear that if 
all propagators involving I2 are pinched, the resulting diagram is the massive double box. 
Explicitly, neglect the second equation in (5.6) and we find that the remaining equation, 

a,2(x)z 2 + a\(x)z + ao(x) = , (5-H) 

which corresponds to the massive double box topology, defines an elliptic curve with the 
genus 1. Then we consider the projection of the variety (5.6) to this elliptic curve, via 

(x,y,z) ^ (x,z) (5.12) 

This is a ramified double covering. The ramified point are determined by the discriminant 
A(z) for the equation, 

b 2 (y)z 2 + bi(y)z + b (y) = , (5.13) 

in y. The simultaneous equations A (2) = and a,2(x)z 2 + a\{x)z + ao(x) = has 8 distinct 
solutions. Then we can use Riemann-Hurwitz formula (2.13) and the fact that the complex 
plane has the genus zero, 

2g(C) - 2 = 2 x (2 x 1 - 2) + 8 • (2 - 1) , (5.14) 

As expected, we get the same conclusion that the curve has the genus 5. Note that 
here we explicitly used the genus information of the massive double-box topology. In this 
sense, Riemann-Hurwitz formula provides an induction relation between the three-loop 
box-penta-box diagram and two-loop double box diagram. 

5.2.1 Degeneracy under specific kinematics 

The degeneracy pattern of the genus-5 Riemann surface is more complicated than that of 
genus-3 Riemann surface of two-loop example. Here we will discuss some degenerate cases. 

If at least one of pi,P2 is massless, Qi(x3,Xi) is factorized and there are two irreducible 
branches. Similarly, if at least one of p^,P5 is massless, Q2(2/3>2/4) has two factors. Also, if 
P3 is massless, Qs(z3, £4) has two factors. When one of pQ,p-?,p8,P9 is absent, the algebraic 
set defined by Qi,i = 1, ..,5 is also reducible, but the factorization is not obvious. Again, 
we can get two irreducible branches from primary decomposition method by Macaulay2 
[37]. 

Let us assume p\ = 0, and the quadratic polynomial of l\ has two factors Qi(xs, X4) = 
/i(x3, £4)72 (£3 j £4)1 where /i,/2 are both linear in £3, £4. The two irreducible branches 
are defined by h = Q 2 , Q3, Q4, Qs) and 7 2 = (f2, Q2, Q3, Q4, Qh)- By the elimination 
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(a) (6) (c) 

Figure 6. Topological pictures of degenerate on-shell equations from planar three- loop box- 
pentagon-box diagram under specific kinematic configurations, (a) the curve is irreducible, and 
the genus is 5, (b) the curve has 2 irreducible branches, each branch has genus 1. They are con- 
nected by 4 points, (c) the curve has 4 irreducible branches, each branch has genus 0. They are 
linked in chain with two adjacent branches connected by 2 points. 

method via Grdbner basis, we can get a plane curve of degree 8 for each branch, so the 
arithmetic genus is qa = 21. There are 8 singular points of the multiplicity \xp = 2, and 2 
singular points of the multiplicity fip = 4, so the geometric genus is go = 21 — 8x1 — 2x6 = 
1. We can also count the intersection points between two curves, and there are 4 single 
points. So the topological picture of this degeneracy can be given by contracting tubes to 
points along dashed lines, as shown in Figure (6.b). 

We consider another kinematic configuration where p\ and p§ are massless. In this 
case both Qi (£3,3:4) and Q2 (2/3, 2/4) are degenerate, and we can get 4 irreducible branches. 
From Macaulay2, we can also get 4 irreducible ideals. For each branch, we can get a plane 
curve of degree 4 via Grobner basis method, so the arithmetic genus is g^ = 3. There are 
3 singular points of the multiplicity fip = 2, so the geometric genus is go = 3 — 3 x 1 = 0. 
Counting intersecting points, we found that there are 8 points in total: the 4 irreducible 
branches are linked in chain, and the adjacent two branches intersect at 2 points. The 
topological picture can be given by contracting tubes to points around dashed lines, as 
shown in Figure (6.c). Each branch is a genus-0 Riemann sphere. 

5.3 Three-loop non-planar box-crossed-pentagon diagram, genus up to 9 

In this subsection, we discuss the three-loop box-crossed-pentagon diagram as shown in 
Figure (7. a). This is a non-planar three-loop diagram with 9 general massive external 
momenta. As usual, the maximal unitarity cut of this diagram gives 11 on-shell equations of 
12 variables, and they define an algebraic curve. Again, there are 3 propagators containing 
only £1, 

D = £j, Dt = (h - Pl f , D 2 = (h - Pl - p 2 f . (5.15) 
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Figure 7. (a) Non-planar three- loop box-crossed-pentagon diagram, (b) Non-planar three- loop 
crossed-crossed-pcntagon diagram. All external momenta are out-going and massive. The loop 
momenta are denoted by £1,^2,^3- 



From on-shell equations Dj = we can parameterize £\ rationally by one free parameter 
x. Similarly, there are 3 propagators containing only £3, 

D = £ 2 3 , D 1 = (£ 3 -p 4 ) 2 , D = (£ 3 - Pl -p 2 -p 3 -p 4 -p 9 ) 2 . (5.16) 

Using on-shell equations Di = 0, £3 is rationally parameterized by one free parameter w. 
However, only 2 propagators 

D = £ 2 , D x = (£ 2 -P6) 2 (5.17) 

containing single loop momentum £2, so £2 is rationally parameterized by 2 free parameters, 
namely, y and z from on-shell equations Di = 0. 

The remaining 3 propagators contain terms of mixed loop momenta, 

£>0 = (*2 " h ~ P5 - P6? , Dl = (h-h-P5-P6-P7) 2 , 

D2 = (£ 1 +£ 3 - Pi -P2- P3 -Pa) 2 ■ (5.18) 

After substituting £i(x), £2(11, z), £^{w) back, they become meromorphic functions. The 
numerators /i,/2,/3 of these 3 meromorphic functions are polynomials in (x,y, z,w), and 
equations 

fi{y,z,w) = f 2 (y,z,w) = h{x,w) = (5.19) 

define the algebraic curve. They are not necessary quadratic. 

Our strategy is to eliminate y and w from the equations and then get a plane curve in 
x and z. This can be done automatically by Grobner basis method. However, it is helpful 
to eliminate y and w step by step, so that we can explicitly show this projection process 
is birational. Furthermore, we can see the induction relation from the two-loop non-planer 
diagram to this diagram. 
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First, note that if we combine the external legs pi, P2, P3 and p$ to one external leg 
and neglect the four cut equations involving l±, then we get the prime case of two-loop 
non-planar crossed-box diagram. In other words, fi(y,z,w) = f2(y,z,w) = defines a 
genus-3 complex curve. Hence we can eliminate y like the case (5.6) and it is birational via 
an inverse transformation like (5.7). The resulting reduced algebraic system is 

agw 8 + aj(z)w 7 + as(z)w 6 + a^,(z)vf' + cn{z)w A + as(z)w 3 + a2(z)w 2 + a\(z)w + ao(z) = , 

b 2 {x)w 2 + b 1 (x)w + b (x) = , (5.20) 

where for generic massive external momenta, Oj, i = 0, . . . 6 and bi, i = 0, . . . 2 are quadratic 
polynomials while clt(z) is linear and as is a constant. Again, we use the resultant to 
eliminate w and get a plane curve F(x, z) = 0. This step is also birational, and the inverse 
map is 

— S H. (Ml) 
<j(x,z) 

where the explicit expressions for p(x, z) and q(x, z) can be found by Grobner basis method. 

Finally, the plane curve F(x, z) = has the degree 20, so the arithmetic genus is 
gA = 171. There are 32 normal singular points of the multiplicity fj>p = 2, one normal 
singular point of the multiplicity fj,p = 4, and one 16- fold point. The 16-fold point has 8 
ordinary tangent lines with the multiplicity 1 and 4 tangent lines with the multiplicity 2. 
Since this 16-fold point is not a normal singular point, we need to perform a blow up. The 
new blow-up curve resolves the 16-fold point into 8 smooth points and 4 normal double 
points. In summary, the genus is 

g G = i(20 - 1)(20 - 2) - 32 • 1 - i • 4 • (4 - 1) - i • 16 • (15 - 1) - 4 = 9 . (5.22) 

So the maximal unitarity cut of the non-planar three-loop box-crossed-pentagon diagram, 
with all massive legs, defines a one-dimensional algebraic curve of the genus 9. 

Alternatively, we can use Riemann-Hurwitz formula to get the genus of the curve 
(5.20). The first equation of (5.20) generates 8 ramified points on the w plane while the 
second equation generates 4 ramified points on the w plane. All these points are distinct, 
so there are 8 • 2 + 4 • 2 = 24 ramified points on the curve (5.20). So Riemann-Hurwitz 
formula reads 

2^- 2 = 4(0 - 2) + 24 , (5.23) 
which again implies that (5.20) has the genus 9. 

5.3.1 Degeneracy under specific kinematics 

Before going to the discussion of degenerate cases, let us first verify an observation of 
on-shell equations 4 . If we forget the 4 equations that containing £±, i.e., Do — D\ = 
Z?2 = D2 = 0, then the remaining 7 equations describe the prime case of two-loop non- 
planar crossed-box diagram, so it should have the genus 3. To see this, we could first 
eliminate 4 variables from 4 linear equations Dq — D\ = Dq — D2 = 0, Dq — D\ = and 



4 We thank Simon Caron-Huot for this observation. 
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I 




Figure 8. Topological pictures of on-shcll equations from the non-planar three- loop box-crosscd- 
pentagon diagram under specific kinematic configurations, (a) the curve has 2 irreducible branches, 
each branch has the genus 3 and they are connected by 4 points, (b) the curve has 2 irreducible 
branches, each branch has the genus 1, and they are connected by 8 single points, (c) the curve has 
4 irreducible branches, each branch has the genus 0, and they are linked in chain, (d) the curve has 
8 irreducible branches, and each branch has the genus 0. 

Dq — Dq — Dq = 0. Then the remaining 3 equations of 4 variables define an equivalent 
curve, and we can birationally project it to a plane curve. The resulting plane curve has 
the degree 8, while there are 12 normal singular points of the multiplicity mp = 2, 1 normal 
singular point of the multiplicity mp = 4, so the genus is 

g G = i(8 - 1)(8 - 2) - 12 • 1 - ±4(4 - 1) = 3 , (5.24) 

which agrees with the observation. 

The genus-9 curve of generic kinematics is complicated. However, as previous examples, 
we can still consider some kinematic limits where on-shell equations are degenerated. In 
these cases the solution space has several branches, and each branch defines a curve with 
the lower genus. Here we will show some simple degenerate cases. 

After solving linear equations, we get 5 quadratic equations, 

Qi(xx,x 2 ) = , Q 2 (yi, 2/2,2i, z 2 ) = , Qz{z ll z 2 ) = , 
Qi{x\,x 2 ,zi,z 2 ) = , Q 5 (yi,y 2 ,zi,z 2 ) = , (5.25) 

where Qi,Q 2 ,Q3 come from and Qa-,Q§ come from the mixed terms £2-^3- 

If pi or p 2 is massless, then Q\ has two factors, while if is massless, then Q3 has two 
factors. In these kinematic configurations, there are two irreducible branches. As usual, we 
can project each branch onto a plane curve. The resulting plane curve has the degree 12, 
so the arithmetic genus is gA = 55. There are 16 normal singular points of the multiplicity 
mp = 2, 1 normal singular point of the multiplicity mp = 4, one 8-fold singular point. 
This 8-fold point is not normal, and after blowing up, we get 2 normal singular points of 
the multiplicity mp = 2. So finally we get, 

g G = 55 - 16 • 1 - ^4(4 - 1) - ^8(8 - 1) - 2 • 1 = 3 . (5.26) 



- 23 - 



Therefore there are two genus-3 curves. We can also calculate the intersecting points 
between these two curves, and find that there are 4 points. The topological picture can be 
reproduced from the genus-9 picture by contracting tubes to points along dashed lines as 
shown in Figure (8. a). 

If pq or P7 is massless, Q2 has two factors. In this case, again, there are two irreducible 
branches. Similarly, after birationally projecting each branch onto a plane curve, we get an 
equation of the degree 12. There are 20 normal singular points of the multiplicity mp = 2, 
one normal singular point of the multiplicity mp = 4 and one normal singular point of the 
multiplicity mp = 8, so finally we get genus 



There are 8 intersection points between the two irreducible branches . So the topological 
picture can be reproduced from the genus-9 picture by contracting tubes to points along 
dashed lines, as shown in Figure (8.b). 

In fact, if we know the genus for generic kinematics, together with the number of 
irreducible branches for specific kinematic configuration and the intersecting points between 
different branches, it is possible to predict the genus of each branch by directly studying the 
topological picture. Take the kinematic configuration in previous paragraph as an example: 
intuitively, these two branches are symmetric, and they should have the same genus. A 
simple calculation determined that there are 8 intersecting points. Then it is clear that 
(8.b) is the only possible topological picture which satisfies these requirements. So we can 
conclude that each branch is genus-1 without doing any further calculation. 

If both pi and pq are massless, Qi and Q2 are factorized, and each of them has two 
factors. Thus, the variety has 4 irreducible branches. The on-shell equation system of 
each branch is quite simple, and the corresponding plane curve for each branch has the 
degree 4. There are 3 normal singular points of the multiplicity mp = 2, so the genus is 
go = (4 — 1)(4 — 2)/2 — 3 = 0. The topological picture is shown in Figure (8.c). 

If we further consider the case with massless pi,P4,Pe, there are 8 irreducible branches. 
The plane curve of each branch has the degree 3, and there is only one normal singular 
point of the multiplicity mp = 2, so the genus go = (3 — 1)(3 — 2)/2 — 1 = 0. After obtaining 
the intersecting points between branches, we draw the topological picture in Figure (8.d). 

5.4 Three-loop non-planar crossed-crossed-pentagon diagram, genus up to 13 

The same methods can be used to study the non-planar three- loop crossed-crossed-pentagon 
diagram, as shown in Figure (7.b), although the computation is more complicated. Nev- 
ertheless we should start from the on-shell equations. There are 2 propagators containing 



9G = 



1(12 - 1)(12 - 2) - 20 • 1 - ^4(4 - 1) - is(8 - 1) = 1 . 



(5.27) 



only t\ 



D =£j, D 1 = {i l -p l f 



(5.28) 



2 propagators containing only I2 




2 



2 



(5.29) 



- 24 - 



and 3 propagators containing only 4 

D =4, D 1 = (e 3 -p 1 -p 2 -p7-p 9 ) 2 , £>2 = (h+P4+P5+P6+P8? ■ (5.30) 

From maximal unitarity cut, we can rationally parameterize 4 by two parameters x, y, 4 
by z,w, and £3 by r. The remaining 4 on-shell equations Do = Di = D2 = D 3 = 0, where 

D = (£ 3 - h - P7 ) 2 , Dx = (4 - h - P7 - P9 f , 

d 2 = (4 + 4 + Pe) 2 , 5 3 = (4 + 4 + Pe + ps) 2 , (5.31) 

become 

fi{x,y,r) = f 2 {x,y,T) = f 3 (z,w,T) = f A (z,w,r) = , (5.32) 

which are polynomial equations. Again, it is clear that this diagram contains two copies 
of non-planar two-loop diagrams. In other words, /1 = /2 = defines a genus-3 curve 
and /3 = /4 = defines another genus-3 curve. Following the previous analysis, we can 
birationally eliminate y and w and get, 

gi(x,r) = , g2{z,r) = , (5.33) 

where g\ is quadratic in x and 52 is quadratic in z. 

We can further eliminate r to get a plane curve F(x, z) = in z. The plane curve 
has the degree 24, so the arithmetic genus is qa = 253. There are 184 normal singular 
points of the multiplicity fip = 2, two normal singular point of the multiplicity \ip = 8. In 
summary, the genus is 

g G = ^(24 - 1)(24 - 2) - 184 • 1 - 1 • 8 • (8 - 1) • 2 = 13 . (5.34) 

Alternatively, each of g± and g2 generates 8 ramified points on the r plane. All these 
points are distinct, so there 8 • 2 + 8 • 2 = 32 ramified points on the curve (5.33). So 
Riemann-Hurwitz formula reads 

2^- 2 = 4(0 - 2) + 32 , (5.35) 

which again implies that the curve from crossed-crossed-pentagon diagram, with all massive 
legs, has the genus 13. 

5.4.1 Degeneracy under specific kinematics 

Again, if we forget the 4 equations containing 4 (or 4)> a s in the box-crossed-pentagon 
example, the remaining 7 equations also describe the generic two-loop non-planar crossed- 
box diagram, thus it should be genus 3. This can be verified by the same method in the 
previous subsection. 

We discuss some kinematic configurations, where the on-shell equations are degener- 
ated. Again, after solving linear equations Dq — D\ = and Dq — Dq — Dq = 0, we can 
formally write the other two quadratic equations as 



Qi(x 1 ,x 2 ,z 1 ,z 2 ) = , Q^(xi,X2,Zx,z 2 ) = 



(5.36) 
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Figure 9. Topological pictures of on-shell equations from the non-planar three- loop crossed-crosscd- 
pentagon diagram under specific kinematic configurations, (a) The curve has 2 irreducible branches, 
and each branch has the genus 3. They are connected by 8 points, (b) The curve has 2 irreducible 
branches, and each branch has the genus 5. They are connected by 4 points, (c) The curve has 
4 irreducible branches, and each branch has the genus 0. They are linked in chain with adjacent 
branches intersecting at 4 points, (d) The curve has 8 irreducible branches, and each branch has 
the genus 0. 



lipi is massless, Qi has two factors Q\ = /1/2, where fi, fi are linear functions. Thus, there 
are two irreducible branches, from ideals I\ = Q2, Q3, Q4, Q5) and I2 = {f2, Q2, Q3, Qi, Q§)- 
These two branches intersect at 8 points, so we can guess that the topological picture should 
be Figure (9. a) and each branch gives a curve with the genus 3. To verify this, as usual, 
we can project the curve onto a plane curve, and the resulting equation has the degree 12. 
There are 40 normal singular points of the multiplicity mp = 2, 2 normal singular points 
of the multiplicity mp = 4, so the genus is 

g G = 1(12 - 1)(12 - 2) - 40 • 1 - 2 • ^4(4 - 1) = 3 , (5.37) 

which is exactly the result shown in Figure (9. a). 

If j>3 is massless, Q3 has two factors. There are two irreducible branches. The corre- 
sponding plane curve of each branch has the degree 12. And there are 38 normal singular 
points of the multiplicity mp = 2, 2 normal singular points of the multiplicity mp = 4, so 
the genus is given by 

g G = 1(12 - 1)(12 - 2) - 38 • 1 - 2 • U(A - 1) = 5 . (5.38) 

The number of intersecting points between two branches is 4, so the topological picture is 
given by contracting tubes to points from the genus-13 picture, as shown in Figure (9.b). 

It is not hard to predict that if both pi,P3 are massless, the topological picture should 
be given by overlapping of Figure (9. a) and (9.b). In other words, there are 4 genus-1 tori, 
which are linked in a chain. In fact, in this kinematic configuration, both Q\ and Q3 are 
factorized, and there are indeed 4 branches. The plane curve of each branch has the degree 
6, so the arithmetic genus is gA = (6 — 1)(6 — 2)/2 = 10. There are 9 normal singular 
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Figure 10. Non-planar three-loop "Mercedes" logo diagram. All external momenta are out-going 
and massive. The loop momenta are denoted by ^i,^2)^3- 

points of the multiplicity rap = 2, so the genus from explicit calculation is 

g G = 10 - 9 ■ 1 = 1 , (5.39) 

as expected. 

If both pi,P5 are massless, Q\, Q2 are factorized. So there are 4 irreducible branches. 
The plane curve of each branch has the degree 6, while there are 10 normal singular points 
of the multiplicity rap = 2, so the genus is 0. If Pi,P3,p$ are massless, we further get 8 
branches. The plane curve of each branch has the degree 3, and there is only one singular 
point of the multiplicity 1, so the genus is again 0. After obtaining the intersecting points 
between branches, we can sketch the topological pictures of both kinematic configurations 
from the genus-13 picture in Figure (9.c) and (9.d). 

5.5 Three-loop Mercedes- logo diagram, genus up to 9 

So far, we only consider the three-loop diagrams with "ladder" type. The same mathemat- 
ical method also works for three-loop diagrams with "Mercedes-logo" topology. 

In this subsection, we explicitly calculate a three-loop "Mercedes-logo" diagram with 
11 propagators and 9 massive external legs as shown in Figure (10). Different from the 
"ladder" diagrams, all possible terms of loop momenta £\£2 , ^1^3, ^2^3 will appear. This 
will somehow complicate the cut equations. There are two propagators containing only £\ 

D = £ 2 , D 1 = (£ 1 -p 1 f , (5.40) 

and three propagators containing only £2 

D = £ 2 , D^^-pe) 2 , D 2 = {£2+Pi? , (5.41) 

and three propagators containing only £3 



D = £ 2 3 , D 1 = (£ 3 -p 3 ) 2 , D 2 = (£3+ Pi ) 2 . (5.42) 



The remaining three propagators contain mixed terms as 



Do = (t 2 -l 1+P7+Pl+P8 f , D 1 = (£ 3 -£ 2 +P4 + P6+Ps) 2 , (5-43) 
D 2 = (£ 1 -£ 3 +p 3 +p 2 ) 2 . (5.44) 

It is not possible to construct linear equations directly from equations of mixed propagators, 
so we can only get 5 linear equations from Di = Di = Di = 0. The one-dimensional curve 
is described by 6 quadratic equations of 7 variables. 

As before, we can parameterize £\ by two variables x and y, £2 by one variable z, and 
£3 by one variables w. The rest equations have the following form 

f{x,y,z)=0 , g(w,z)=0 , h(x,y,w) = . (5.45) 

We birationally project the variety on a plane curve F(y, w) = 0. This plane curve 
has the degree 20, 36 normal double point, one normal 16-fold point and one normal 4-fold 
point. So the geometric genus is, 

(20-l)(20-2) 16(16 - 1) 4(4 - 1) 
9G = 2 36 2 2 = ' ^ ' 

Interestingly, genus of Mercedes logo diagrams equals to the genus of box-crossed- 
pentagon diagrams, they both have genus 9. We examine some simple degenerate equa- 
tions under specific kinematic configurations, and find that the topological pictures are 
the same as box-crossed-pentagon. If P4 is massless, there are two irreducible branches. 
After birationally projecting each branch onto plane curve, we get an equation of degree 
12. There are 18 normal singular points of mp = 2, 1 normal singular point of mp = 4 
and 1 normal singular point of mp = 8, so the genus is 3. These two branches intersect 
at 4 points, and the topological picture is the same as shown in Figure (8. a). If p\ is 
massless, cut equations also degenerate to two irreducible branches. The plane curve of 
each branch has degree 12. There are 20 normal singular points of mp = 2, 1 singular 
point of mp = 4 and 1 singular point of mp = 8, so the genus is 1. There are 8 intersecting 
points between two branches, and the picture is the same as Figure (8.b). If both pi,p$ are 
massless, there are 4 irreducible branches. The plane curve after birationally projecting 
each branch has degree 6 and 10 normal singular points, so each of them is genus-0 curve. 
The topological picture is the same as Figure (8.c). Further considering pi,P4,Pe massless, 
we get 8 irreducible branches. Each branch can be projected onto a plane curve of degree 3 
and 1 singular point, so they are again genus-0 curve. After getting the intersecting points 
among branches, the topological picture can be found to be the same as the one shown in 
Figure (8.d). 



6 Conclusion 

In this paper, we study the global structure of on-shell equations from generalized unitarity 
cut of loop diagrams. We focus on the L-loop diagram with 4L — 1 unitarity cuts in four- 
dimensional theories. There is one degree of freedom for the solution space, which is a 
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complex algebraic curve or union of complex algebraic curves. Since the topology of a 
complex algebraic curve is completely determined by the genus, in order to study the 
global structure, we compute the genera of algebraic curves. This is systematically done 
by computational algebraic geometry methods. 

In this paper, we used two algebraic geometry methods to determine the genera of 
curves from unitarity cuts: (1) the relation between arithmetic genus and geometric genus 
(2) Riemann-Hurwitz formula. These methods clearly work for all algebraic curves, and 
are quite efficient for those from two and three-loop unitarity cuts. 

The only algebraic curve for the one-loop case is from the one-loop triangle diagram 
with 3 propagators. The defining equation is simply a conic section, and it is equivalent 
to a genus-0 Riemann sphere, so we can always find the rational parametrization for the 
solution space. 

Two-loop diagrams with 7 propagators also define algebraic curves. We find that the 
genus could be as high as 3, in the two-loop non-planar crossed-box diagram with the 
maximal number of massive external momenta. In the degenerate limit, i.e., some external 
momenta being massless or absent, the solution space degenerates to the union of several 
complex curves. These curves have lower genera, 1 or 0. As the global structure analysis 
in [23], from intersection structure of these curves, we can show that the union of these 
curves exactly comes from a degenerate genus-3 curve. 

The topology of algebraic curves from three-loop diagrams with 11 propagators is more 
complicated. We consider four examples. The pentagon-box-box diagram is, in some sense, 
similar to the two-loop double box diagram. The curve from the planar box-pentagon-box 
diagram with generic kinematics has the genus 5. For the non-planar box-crossed-pentagon 
diagram with generic kinematics, the curve can have the genus 9. And for the non-planar 
crossed-crossed-pentagon diagrams, the curve can have the genus 13. Again, for diagrams 
with massless external legs, we show that the global structure of solutions comes from a 
degenerate high-genus curve. 

The information of genus and the topological picture under degenerate limit is impor- 
tant for calculating multi-loop amplitude via unitarity cut method: 

• The genus is the criterion for rational parametrization: if the curve has the genus 
zero, then it can be rationally parameterized. Otherwise, the rational parametrization 
does not exist. This property is important, since for unitarity computation, we may 
parameterize solution space before tree amplitude computations. So it is useful to 
know the difficulty of parametrization before unitarity computation. 

• The topological picture under degenerate limit is important for calculating expansion 
coefficients of integrand basis. Usually, all branches should be studied to get these 
coefficients. Knowing the genus for generic kinematics, it is possible to predict the 
topology for branches under degenerate limits from the knowledge of intersection 
pattern. Especially, we can use the genus information to check if all the cut solutions 
are identified or not. 
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We comment that in all cases considered in this paper, the genera of curves, if not 
zero, are always odd. We expect that this feature can be understood by the loop-by-loop 
induction relation, in the future work. 

Unitarity cuts with fewer propagators are more complicated, since in these cases, the 
on-shell equations define algebraic (hyper-)surfaces instead of curves. The topological 
structures of surfaces are much more complicated than those of curves. However, we 
expect that the techniques of computational algebraic geometry would still be important 
for obtaining the topological information of unitarity cuts. Analysis of the global structure 
of generalized unitarity cuts is just the beginning of multi-loop amplitude calculation, and 
more works, such as parametrization or the branch-by-branch polynomial fitting [34], can 
be done based by this information. 
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